Numerical Simulation of an Electroweak Oscillon 
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Numerical simulations of the bosonic sector of the SU{2) x U(l) electroweak Standard Model in 
3 + 1 dimensions have demonstrated the existence of an oscillon — an extremely long-lived, localized, 
oscillatory solution to the equations of motion — when the Higgs mass is equal to twice the W ± 
boson mass. It contains total energy roughly 30 TeV localized in a region of radius 0.05 fm. A 
detailed description of these numerical results is presented. 

PACS numbers: 11.27.+d 11.15.Ha 12.15.-y 

INTRODUCTION 



While static, localized soliton solutions to the equations of motion of nonlinear field theories have been well studied, 
and are of interest in many applications 0, 0| , no known examples exist in the electroweak Standard Model (although 
>— * ' there do exist extended electroweak string solutions HQ). However, much less is known about the existence of 
localized solutions that oscillate in time, known as breathers or oscillons. (The latter term was originally introduced 
' to describe similar phenomena in plasma physics Q.) In some models, such as the sine-Gordon breather [y] and Q-ball 
0, one can use conserved charges to prove the existence of exact, periodic solutions. But oscillons have also been 
found in many nonlinear field theories that do not contain either static solitons or conserved charges. These solutions 
either live indefinitely or for extremely long times compared to the natural timescales of the system, 
i' For scalar theories in one space dimension, oscillons have been found to remain periodic to all orders in a perturbative 
expansion [fj and are never seen to decay in numerical simulations [8|, but can decay after extremely long times via 
i ' nonperturbative effects Q or by coupling to an expanding background [Toj ] . In both </> 4 theory in two dimensions 
11 , [H | and the abelian Higgs model in one dimension [l3| and in two dimensions [3] , oscillons have been found that are 



not observed to decay. In </> 4 theory in three dimensions, however, one finds long-lived quasi-periodic solutions whose 
lifetime depends sensitively on the initial conditions 15, 16, 17, lH, Similar behavior is present in other scalar 



theories in three dimensio ns Ipoll a nd in higher dimensions [2 If . Phenomenologically, small Q-balls were considered as 
_— , , dark matter candidates in |22j, |23j, |24|, |25j , axion oscillons were considered in | 2 6| . and the effects of oscillons and other 
aspects of nonequilibrium dynamics in and after inflation were studied in [27U281. l29|. Oscillons and related solutions 
have also been studied in connection with phase transitions [1(3] , monopole systems [3lj , QCD [32| , and gravitational 
<^> ' systems [Hj|. 

C ■ . Recent work 34{ demonstrated numerically the existence of an oscillon in the bosonic sector of the electroweak 
Standard Model, when the mass of the fundamental Higgs is exactly twice that of the gauge bosons. (A similar 
mass relation also arises in the study of embedded defects [35j.) This result was based on previous work (3^|, which 



found oscillons in spontaneously broken pure SU(2) Higgs-gauge theory with the same 2 : 1 mass ratio. In that 



model, one can consider field configurations restricted to the spherical ansatz 37J, meaning they are assumed to be 
invariant under combined rotations in space and isospin, also known as grand spin rotations. Within this ansatz, 
the system can be described by an effective theory of fields depending only on r and t, which greatly simplifies the 
numerical analysis. In [iiij this numerical simulation was extended to a fully three-dimensional spatial lattice with no 
assumptions of rotational symmetry, making it possible to also include the U(l) hypercharge field (which breaks the 
grand spin invariance of the spherical ansatz). The resulting simulation comprises the full electroweak sector of the 
Standard Model without fermions. Here we extend that analysis and describe its results in more detail. We use the 
same SU (2) gauge coupling g and Higgs self-coupling A as in the pure SU (2) theory, meaning that the Higgs mass is 
twice the mass of the bosons, and set the U(l) coupling g' so that the mass of the Z° boson matches its observed 
value. 

Ong oin g an alytic work [38| has shed some light on the 2 : 1 mass ratio by using a small amplitude approximation 
fill. Ilfj, 32. 3t| to construct oscillons in a simplified version of the spherical ansatz theory. In this analysis, one begins 



by assuming that each field in the oscillon profile has large width, so that at large distances it falls like exp(— emr) 
where m is its mass. There, the amplitude is small and the oscillations obey a linear dispersion relation, which implies 
u> = mvl — e 2 . The linear, dispersive gradient terms in the equation of motion are then of order e 2 . They must 
be balanced by nonlinear terms to obtain a stable solution. Since the leading nonlinearity is typically given by a 
quadratic term in the equations of motion, this requirement implies that the field amplitudes must be proportional 
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to e. In a multiple-field model, one must also ensure that the terms giving interactions between different fields are 
resonant with the dispersive linear terms, so that their effects are not washed out over many cycles. As shown in 38], 
the 2 : 1 mass ratio arises naturally in this analysis: Since the fields' oscillation frequencies are tied to their masses, 
imposing a resonance condition on their frequencies is equivalent to fixing a particular mass ratio. Although this 
analysis has so far only been carried out in simplified models, we will see below that the oscillon observed numerically 
in the full electroweak theory is of small amplitude and large width, so similar techniques are potentially applicable 
in this case as well. 

In all known oscillons, each field oscillates with a frequency below its mass, so that it couples to dispersive linear 
waves (which have lu = \/k 2 + m 2 > m) only through nonlinear interactions. The fields then converge to a configu- 
ration in which this decay channel is also suppressed. Because the electroweak theory includes the massless photon 
field, which can radiate in arbitrarily low frequencies, one might expect the oscillon to decay rapidly by emitting 
electromagnetic radiation, but it does not. Instead, after initially shedding some energy in this way, the system settles 
into a localized solution that no longer radiates and remains stable for as long as we can follow it in numerical sim- 
ulations. In preliminary work that provided motivation for the current investigation, similar behavior was observed 
both when an additional massless scalar field was coupled to oscillons in onc-dimensional (f> A theory and when an 
additional spherically symmetric massless scalar field was coupled to oscillons in the spherical ansatz model. In each 
case, after shedding some energy into the massless field, the oscillon arranges itself in a neutral configuration that 
no longer couples to the massless field. This mechanism may be similar to the suppression of nonlinear coupling to 
dispersive waves that is common to all oscillons. 



CONTINUUM THEORY 



We begin from ST/ (2) x U{\) electroweak theory in the continuum, ignoring fermions, and follow the conventions 
of (40| . The Lagrangian density is 



C = --F^F» V - -F^ v ■ + {D^D»<f> - A(|$| 2 - v 2 ) 2 



(1) 



where the boldface vector notation refers to isovectors. Here $ is the Higgs field, a Lorentz scalar carrying U(l) 

hypercharge 1/2 and transforming under the fundamental representation of SU(2). The metric signature is H . 

The SU{2) and 17(1) field strengths are 



Ffj.v = d^Wv - d v W» - gW^ x W v , = 3^B V - d v B^ , 

and the covariant derivatives are given by 



= I d„ + i 9 -B^ + i|r • W^J $ , D^F^ = d»F^ - gW» x F^ . 

where t represents the weak isospin Pauli matrices. We obtain the equations of motion 

= J v , D^F^ = J v , D"D^ = 2X(v 2 - |$| 2 )$ , 

where the gauge currents are 

J v = g'Im(D 1/ $y$, J u =gIm(D 1/ <f>y T $. 
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(2) 



(3) 



(4) 



(5) 



We work in the gauge Bq = 0, Wo = 0. With this choice, the covariant time derivatives become ordinary derivatives 
and we can apply a Hamiltonian formalism. The energy density is 



J=x,y,z 



B) + W j IV ^ (F 2 j + F kj ■ F kj ) 
k>j 



\M 2 + £ (D^(D^) + X(\^\ 2 



(6) 



]=x,y,z 



where dot indicates time derivative. The integral over space of this quantity is conserved by the time evolution. From 
the equations for Bq and Wq, we obtain the Gauss's Law constraints, 



£ ^4-70 = 0, J2 0iWWo = O, 



(7) 



J=x,V,z 
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where the charge densities are 

J =.g'Im$ t $, J =5lm4 t r$. (8) 

These constraints remain true at all times, at all points in space, assuming they are obeyed by the initial value data. 

Although the numerical calculation will be done using the underlying gauge fields W M and B^, because of spon- 
taneous symmetry breaking the physical content of the theory is better described by the fields of definite mass and 
electric charge 

W± = ±=[{W»-x)±i(W»-y)] , 
Zl = (W M - z)casO w -B ll smOw, 

An = cos e w + (Wn ■ z) sin 6 W , (9) 

where x, y, and z denote unit vectors in isospin space and 9w = arctan^'/fl 1 ) is the weak mixing angle. The 
fields have mass mw = gv/y2 and electric charge ±e = ±g'cos8w, the Z® field has mass mz = mw/ cosOw and 
zero electric charge, and the photon field has zero mass and zero electric charge. The only other physical degree 
of freedom in the theory is the magnitude of the Higgs field, with mass mu = 2vV~\ and zero electric charge. 



LATTICE THEORY 



To analyze the classical equations of motion numerically, we use the standard Wilsonian approach [4l| for lattice 
gauge fields (for a review see j42|), adapted to Minkowski space evolution as in 43, 44, 45| . The U(l) and SU(2) 
gauge fields live on the links of the lattice and the Higgs field lives at the lattice sites. We use a regular lattice with 
spacing Ax and determine the values of the fields at time t+ = t + At based on their values at times t and t—=t— At. 



Throughout, we will use the same notation and conventions as 34 1 
We associate the Wilson line 



JJP _ e ig'B 1 ' ] &x/2 e igWj- T Ax/2 



(10) 



with the link emanating from lattice site p in the positive j th direction. We define the Wilson line for the link emanating 
from lattice site p in the negative j direction to be the adjoint of the corresponding Wilson line emanating in the 
positive direction from the neighboring site, U^j — {U^~ ] )\ where the notation p± j indicates the adjacent lattice 
site to p, displaced from p in direction ±j. At the edges of the lattice we use periodic boundary conditions. 
The equation of motion for the Higgs field at site p is 



$ p (i+) = 2<P p {t) - $ p (£_) + At 2 <t> p (t) , 



(11) 



where 



E 

j = ±x,±y,±z 



U p (t)<i>P +: >(t) - &>(t) 

Ax 2 



For the gauge fields, we have 



2A (v 2 - \<S> p {t)\ 2 ) <$> p (t). 



(12) 



Uf(t+) = exp 



^(*),(13) 



where 



U W)W = u?wp J W^ (t)u^ (t) 



and 



Jj = 9 Im 



, $f(t)tC/?(t)$f+J(t) 



<$>p(t)Wm{t)<f>p+j(t) 



Jj = g Im ■ 

Ax J Ax 



are the gauge currents. Here we have defined the logarithm of a 2 x 2 matrix in the form of Eq. (fT0|) as 

logUV = ^(g'B? + gW?-T), 



(14) 



(15) 
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which gives the more familiar gauge fields in terms of the link variables. We note that \ogXY ^ logX + log Y" when 
the matrices do not commute. 

The U(l) and SU(2) matrices in Eq. (fTTTj) are stored separately in the numerical code. To represent the f7 (1) matrix 
Ui = e lB , just the real quantity 9 = g'B p Ax/2 is actually stored. Any SU(2) matrix can be written as 



Xi x 2 



(16) 



so only the two complex elements of the top row need to be stored. (This representation is redundant, since \x\\ 2 + 
\x 2 \ 2 — 1, but more efficient computationally than storing three real quantities and reconstructing the fourth.) The 
logarithms and exponentials needed to convert between the group and the algebra can be computed efficiently using 



in - t sin ( 



cos 9 + in z sin 9 in x sin 9 + n y sh 
in x sin 9 — n v sin 9 cos 9 — ih z sin ( 



(17) 



where n is a unit vector and the link matrices have fi9 = W p gAx/2. 

We note that this discretization differs slightly from the standard approach used in 43|, 44, 45|. In our language, 



their discretization is equivalent to replacing sin # — > # and cos 9 — ► \/l — 9 2 when computing both the logarithm and 
the corresponding exponential. While the approach we are using corresponds a little more directly to the continuum 
equations, any differences are of higher order in the lattice spacing. Numerical experiments show that their approach 
yields completely equivalent results, and is somewhat more efficient computationally, since it avoids the need to 
compute trigonometric functions in this conversion. 
The energy density at p is then 



|exp (log Uf(t+)- log Uf(t-)) | 



(2Ai) : 



£ 

i'>3 



Ax 2 



(2At) 2 



whose integral over the whole lattice is conserved. Here we have defined 



ll2 TrlogC/; 

\lJ p \\ = 3 
I ill g a Ax 2 



TrrlogC/J) 1 • (TrrlogC/J 
g 2 Ax 2 



= \B?\ 2 



for any U{2) link matrix. 

At every lattice point, Gauss's Law, 



E 



iog^(t + )^(t)t + io g £/; 3 .(t + )^ J .(t)t 

2iAx 2 At 



0. 



is also maintained throughout the evolution, where the charge densities are given by 

( <$> p {t+) - <f> p {t)\ f 



wM' "H: t ' (t) V«'M 



J = glm 



At 



T$ P (t). 



(18) 



(19) 



(20) 



(21) 



This requirement will provide a stringent check on the correctness of the numerical simulation. Here we have computed 
Gauss's Law at time t + At/2, which is obeyed exactly by the discrete equations of motion for any time step and 
lattice spacing. In [34j ]. Gauss's Law at time t was used; it is only obeyed to order At 2 , but as a result it also provides 
a rough estimate of whether the time step is small enough. 



SPHERICAL ANSATZ 

With the U(l) field included, the grand spin symmetry of the spherical ansatz used in [3(| is broken and field 
configurations will not maintain this symmetry under time evolution. The continuum theory does still preserve 
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invariance under grand spin rotations around the z-axis, but the Cartesian lattice provides a small breaking of all 
rotational symmetries. As a result, field configurations that start within the spherical ansatz are not constrained to lie 
in any reduced ansatz at later times. (We will also demonstrate the oscillon's stability under explicitly nonspherical 
deformations below.) Nonetheless, because we will use the spherical ansatz as a starting point to obtain our initial 
conditions, it will be helpful to analyze it in more detail. We will see that the electroweak oscillon retains much of 
the structure it inherits from these initial conditions. 

For our choice of gauge, the spherical ansatz takes the form [37| 



W S = 



1 

9 
1 

9 



a(r, t) 



ai(r, t)r ■ rf^ 
\p(r, t) — iv(r, t) t ■ r 




1 



\ 7(7, t),~ \ 
rr,) — (r x t) 3 



(22) 

where r is the position vector, r = \r\ is the distance from the origin, and f = r/r is the unit radial vector. 
Configurations in this ansatz are then described by reduced fields a\, a, 7, /x, and ^, all of which depend only on r 
and t. The field definitions have been chosen so that the reduced fields match those used in (36|, even though the 
conventions for the three-dimensional theory used here are slightly different. 

These configurations are in the grand spin zero channel, meaning they are symmetric under simultaneous rotations 
in space and isospin. The gauge field Wj has isospin i = 1 and internal angular momentum s = 1. These two spins 
can be coupled together to yield total generalized angular momentum 0, 1, and 2. To obtain grand spin G = 0, 
these combinations must then be coupled with equal orbital angular momenta I = 0, £ = 1, and 1 — 1 respectively, 
corresponding to monopole, dipole and quadrupole spatial distributions. These three possibilities are reflected in Eq. 
(|2"2")l through the three terms a(r, t), 7(7-, t), and ai(r, f). 

We have written the Higgs field as a matrix times a fixed isospinor. This matrix transforms under both the gauged 
SU(2) L and global SU(2) R isospin transformations. (We are only considering global rotations in both cases, however.) 
Under both transformations it has isospin i = 1/2, giving total isospin i = or i = 1. Since the Higgs is a Lorcntz 
scalar, with zero internal angular momentum, to obtain G — these two possibilities must be coupled to £ — and 
I = 1 respectively, corresponding to monopole and dipole spatial distributions. These possibilities appear in Eq. (|2"2"]) 
as the terms n(r,t) and v(r,t). 

Although the spherical ansatz does not contain the U(l) field, to leading order in 0w we can find the electric charge 
density created by a spherical ansatz configuration for our choice of gauge [4t| . 

2ez 

Jo = ~3-~2 (7" - a l) ■ ( 23 ) 

r°g 

The charge shows a dipole structure centered on the z axis — as we would expected since the electromagnetic 
interactions break the grand spin symmetry by selecting the z direction in isospin. We note that this electric charge 
density is time independent (and thus does not radiate) if the a and 7 fields vary sinusoidally in time with the same 
frequency. 

NUMERICAL SIMULATION 



The initial conditions for the simulation are obtained starting from an approximate functional fit to the solutions 
that were found in S'[/(2)-Higgs theory using the spherical ansatz (3^. These results, with slight modifications, 
provide the initial data for the Wj and $ fields, and the initial Bj field is chosen to vanish. In order to guarantee that 
the initial configuration obeys Gauss's Law in the full SU(2) x U(l) theory, we generate the spherical ansatz fit at a 
point in the cycle where the time derivatives are smallest, and then set all time derivatives to zero. We note that in 
pure SU (2) Higgs-gauge theory, this restriction would not be necessary, because even though an approximate fit with 
nonvanishing time derivatives will not obey Gauss's Law, we can restore Gauss's Law by adjusting $(£+) slightly via 
an SU(2) transformation at each point, 



$ncw(i+) = 



$old(*H 



$(t) 



U p <S>(t) , 



with 



U p = exp 



E 



log U[ (t, 



)^ft)t + io g t/rft + )rjr(t)t 



ff 2 Az 2 |$ old (;+) ||<&(t0 1/2 



(24) 



(25) 
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This procedure has been used successfully to reproduce spherical ansatz solutions with nonvanishing time derivatives 
at t = in a fully three-dimensional simulation of pure SU{2) Higgs-gauge theory, but it cannot be extended to the 
SU(2) x U(l) theory because carries both charges, and thus cannot be adjusted to satisfy both constraints at once. 
Therefore we will consider only initial conditions in which all fields have zero time derivatives, so that Gauss's Law is 
trivially satisfied. 

To construct the initial conditions, we begin from the spherical ansatz form of Eq. (|22p . We work in units where 
v = l/\/2. Since we are dealing with purely classical dynamics, we can rescale the fields to fix the SU(2) coupling 
constant at g = y2, so that the W ± mass is then mw = gvj\[2 = l/v2. With this rescaling, we must also introduce 
an overall factor of g 1 1 g%j multiplying the total energy, where gw — 0.634 is the true weak coupling constant. (This 
factor was incorrectly omitted in the original version of (23|.) We choose A = 1, so that the Higgs mass is twice the 
W mass, run = 2vV\ = V2. Finally, we fix g' = 0.773, so that the ratio g'/g matches its observed value and the Z° 
boson has the correct mass. With these choices, one unit of energy is 114 GeV, one unit of time is 5.79 x 10 -27 sec, 
and one unit of length is 1.74 x 10~ 18 m. In these units, we take the following initial configuration for the radial fields, 

ai (r) = x(0.H7x + 0.016xr)(sech2xr) 1/8 , 

fi(r) = 1 - 0.138xsech-^- , 
6.75 

yr 

v{r) = 0.017xfsech — , 
5 

a(r) = 0.117x 2 rsech^ , 

7 (r) = 0, N (26) 

where the adjustable parameter x allows us to include a combined rescaling of the fields' amplitudes and r-dependence, 
as is commonly used in a small amplitude analysis @, I32I . While x = 1 gives an approximation to the spherical 
ansatz solution of [36}, a slightly larger value appears to be necessary for the configuration to settle into a stable 
solution in the full SU{2) x U(l) model. Here we will use x = 1-15. The first term in parentheses in the definition 
of ai(r) is scaled with an additional x so that it matches the coefficient of a, ensuring that a, a\ — a/r, 7/r, and 
v all vanish as r — > 0, as required for regularity of the fields at the origin. Within the spherical ansatz simulation, 
these initial conditions converge to a long-lived oscillon in the pure 5C/(2)-Higgs theory, which is never observed to 
decay. As a check of the numerical calculation, the full three-dimensional simulation agrees with the spherical ansatz 
simulation when the U(l) interaction is turned off. 

Although initial conditions of this form do settle into stable oscillon configurations in the SU(2) x U(l) theory, it 
is helpful to make a minor modification to them that is outside the spherical ansatz: setting the r z -component of Wj 
to zero brings the initial conditions significantly closer to the localized solution that the fields ultimately converge to. 
While we obtain an equivalent oscillon solution in both cases, this modification reduces the energy shed as the oscillon 
forms. Doing so provides a significant technical benefit, because the radiation emitted as the configuration settles 
into the oscillon solution can wrap around the periodic boundary conditions, return to the region of the oscillon, and 
potentially destabilize it. To avoid this problem, the energy density in this radiation, which spreads throughout the 
volume of the simulation, must be small compared to the oscillon's energy density. As long as the lattice volume is 
large enough compared to the oscillon size, this radiation is sufficiently diffuse that it does not affect the oscillon's 
evolution. We use a lattice of size L — 144 on a side in natural units, which is more than enough to satisfy this 
criterion. For L > 100, changing the lattice size simply changes the pattern of noise caused by electromagnetic 
radiation superimposed on the oscillon region, but does not affect oscillon properties or stability. We can therefore be 
certain that there is no coherent structure to this unphysical radiation that could possibly be necessary for the oscillon's 
stability. Its only potential effect is to destabilize the oscillon, and it only does so when artificially concentrated by a 
small lattice (e.g. of size L < 100). In numerical experiments, these destabilization effects are actually much weaker in 
the electroweak model than in pure scalar or SU(2) Higgs-gauge models, because in the electroweak model the radiated 
energy ends up almost entirely in the electromagnetic field, while the oscillon arranges itself to be electrically neutral. 
For this reason, it is not necessary to use absorptive techniques such as adiabatic damping [ll| or an expanding 
background although both have been applied successfully to this problem as well. However, clearly it is helpful 
to adjust the initial conditions to be as close as possible to the true oscillon configuration, to minimize the amount 
of unwanted energy emitted as the configuration settles into the oscillon solution, and therefore limit the numerical 
costs associated with a larger lattice. 

Starting from the modified spherical ansatz initial conditions, we let the system evolve for as long as is practical 
numerically, and see no sign of oscillon decay. We use lattice spacing Ax — 0.75, though Ax = 0.625 and Ax = 0.25 
were verified to give completely equivalent results in correspondingly smaller tests. The time step is At = 0.1. Time 
steps of 0.05 and 0.025 also gave equivalent results, although in this case one must take into account the fact that 
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FIG. 1: Energy in a spherical box of radius 28 as a function of time in natural units. The initial conditions are given by the 
modified spherical ansatz form given in the text, in which the t z component of the gauge field is set to zero, with \ = 1-15. 
Two values of the Higgs self-coupling A are shown. For A = 1, the masses of the Higgs and W fields are in the 2 : 1 ratio needed 
for oscillon formation and the solution remains localized throughout the simulation. Here one unit of energy is 114 GeV, one 
unit of time is 5.79 x 10~ 27 sec, and one unit of length is 1.74 x 10 -18 m, giving a total energy of roughly 30 TeV within the 
box radius of roughly 0.05 fm. A transient beat pattern is also visible. For A = 0.95, the mass ratio is 1.95 : 1. In that case, 
there is no stable object and the energy quickly disperses. 




!,500 



time time 



FIG. 2: Energy in the spherical box for a variety of values of A. For A = 1, the Higgs mass is twice the W ± mass and no decay 
is observed. When the Higgs mass is just below this value, we see a region of meta-stability. For A < 1, the fields decay by first 
collapsing inward before dispersing, while for A > 1 the fields simply disperse outward. 



this change also slightly alters the initial conditions: To set the initial time derivatives to zero, the simulation sets 
the first two time slices equal. Changing the time step thus changes the time at which the field configuration matches 
its value at t = 0, representing a slight perturbation of the initial conditions. This change slightly alters the initial 
transient behavior as the fields approach the oscillon, but these differences quickly disappear and the simulations 
approach equivalent oscillon configurations. 

Total energy is conserved to a few parts in 10 3 for At = 0.1, which improves with At 2 as expected for our second- 
order algorithm. We check Gauss's Law by monitoring the left-hand side of Eq. ([20| . which we verify vanishes to 
machine precision throughout the simulation. 1 It is necessary, however, to use double precision to avoid gradual 



1 One can instead evaluate Gauss's Law at time t instead of t + At/2 as in |'i4l . In that case, we square the left-hand side of Eq. (120 ft . 
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FIG. 3: Decay of the oscillon for A / 1. One of the gauge fields and the magnitude of the Higgs field at the origin are shown 
as functions of time. In the left panel A = 0.99375 and the oscillon decays by collapsing inwards, creating a large amplitude 
fluctuation at the origin before dispersing. In the right, panel, A = 1.00625, and the oscillon decays by expanding outwards. 



degradation in this result. For the parameters as given above, a run to time 10, 000 takes roughly 40 hours using 24 
parallel processes, each running on a 2 GHz Opteron processor core. 2 
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FIG. 4: Energy in the spherical box for two different box radii in the simulation of Fig. [T] with A = 1. The transient beat 
pattern represents a "breathing" perturbation in which the oscillon stretches and compresses slightly. For the larger box size, 
less energy flows in and out of the box during this process, and so the observed beat amplitude is smaller. 

Fig. [1] shows the energy in a spherical box of radius 28 as the fields are evolved from these initial conditions. When 
the Higgs mass is twice the mass, a small amount of energy is initially emitted from the central region, with 
the rest remaining localized for the length of the simulation. If the masses are not in this ratio, however, the initial 
configuration quickly disperses. Fig. [2] shows the growth in oscillon lifetime as A approaches this critical value. We 
see a region of meta-stability when the Higgs mass is just below the 2 : 1 ratio. For A < 1, the fields first collapse 
toward the origin before dispersing, while for A > 1 they simply spread outward. This behavior is shown in Fig. [3] 
Other "special" ratios, such as mjj = 2mz, did not form stable objects from these initial conditions. 



take its trace, and then take the square root of the result. For a typical run with At = 0.1, the integral of this quantity over the lattice 
never exceeds 0.025 and shows no upward trend over time. For smaller At, we see the expected 0(At 2 ) improvement in this result. 
2 The parallel CH — h code used for these simulations is available from http://community.middlebury.edu/~ngraham 
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The spherical box contains approximately 3% of the total volume available to the simulation. Its radius has been 
chosen to be just large enough to enclose nearly all of energy density associated with the stable oscillon. As a result 
of this choice, the A = 1 graph also shows a transient beat pattern. It represents a "breathing" or "ringing" motion, 
in which the oscillon gradually expands and contracts slightly over many periods, accompanied by a corresponding 
modulation of the field amplitudes. This process causes a small amount of the oscillon's energy to move in and out 
of the box. As we would expect, when a larger box size is used, the "breathing" is more completely contained within 
the box and the graph of the energy in the box flattens out, as shown in Fig. [4] Similar beats appear in the SU (2) 
spherical ansatz oscillon [36j | , but in the electroweak oscillon their amplitude decays much more rapidly. 
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FIG. 5: A snapshot of the gauge fields in the x — plane for the simulation of Fig. [T]at time t = 50, 000. Subscripts refer to 
spatial components. 

To illustrate the field configurations that make up the oscillon, we graph the fields at time t = 50, 000 for the 
two-dimensional slice x — 0. Fig. [5] shows the gauge field components. It is most illustrative to consider a linear 
superposition of the fields, as shown in the figure. Fig. [5] shows the electric fields, which are given by the time 
derivatives of the gauge fields for our choice of gauge. Fig. [7] shows the components of the Higgs field and its first time 
derivative, and Fig. [8] shows the magnitude of the Higgs field and the first time derivative of this quantity, together 
with the total energy density. The oscillon is constructed primarily out of the lower component of the Higgs field, the 
imaginary part of the upper component of the Higgs field, the x and y spatial components of the Wj fields, and the z 
spatial component of the Z® field. We see the multipole structures we anticipated from the spherical ansatz analysis. 
The Higgs field contains monopole and dipole fluctuations. The photon field Aj contains delocalized background 
radiation that was emitted as the oscillon formed from the initial conditions. As we would expect from Eq. (|23[) . it 
has a dipole structure. In the spherical ansatz, the Wj and Z® fields can potentially contain monopole, dipole, and 
quadrupole components. Here we see significant monopole and quadrupole structures, but only a very small dipole 
component, which appears in Z®. As a result, the electric charge we estimate from Eq. (|23[) is very small, as is the 
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FIG. 6: A snapshot of the electric fields (time derivatives of the gauge potentials) in the x — plane for the simulation of Fig 
Q]at time t = 50, 000. Subscripts refer to spatial components. 
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FIG. 7: A snapshot of the Higgs field and its time derivatives in the x = plane for the simulation of Fig.[T]at time t = 50, 000. 
Subscripts refer to components of the Higgs field. 



true value from the numerical simulation; the oscillon is decoupled from the electromagnetic background. 3 

Each excited field oscillates at a frequency just below its mass. In our units, these oscillations have typical amplitude 
of order 0.1 and typical radius of order 10. By comparing the total number of cycles to the total time, we find 



3 While the multipole analysis is instructive as a description of the field configuration, is is important to note that because the oscillon 
has large spatial extent compared to its period of oscillation, it is in exactly the domain where the standard multipole expansion for the 
electromagnetic radiation emitted is invalid. 
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FIG. 8: Left panel: A snapshot of the magnitude of <j> an d its first time derivative in the x — plane for the simulation of 
Fig. [T]at time t = 50, 000. Right panel: A snapshot of the energy density the x — plane for the simulation of Fig. [T]at time 
i = 50,000. 
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FIG. 9: Oscillon fields at the origin as functions of time. The left side shows one component of the SU(2) gauge field. The 
upper graph shows the full extent of the simulation. On this scale, the individual oscillations are too small to be seen. Instead, 
we see the decaying beat pattern from the transient "breathing" motion. The lower graph shows the oscillation of the field for 
a short time at the end of the simulation (when the transient effects have decayed away) . The right side shows the magnitude 
of $ in the same way. It oscillates with fundamental frequency twice that of the gauge field. 

(jji = 1-404 for the Higgs field components and low = 0.702 for the gauge field components. These properties are all 
very similar to the spherical ansatz oscillon. They are also consistent with a small-amplitude analysis, as described 
in the Introduction, with e of order 0.1. In Fig. [51 oscillon fields at the origin are shown as functions of time. The 
fundamental oscillation of each field is modulated by the decaying beat pattern. 

The oscillon we have seen is not significantly altered by small perturbations of the initial conditions. As an 
example, in Fig. [10] we show the results of a run in which the rotational symmetry has been explicitly broken. We 
take initial conditions as before, except we introduce different rescalings of the x, y and z coordinates in the definition 
of r. As an additional numerical check, this run also uses a smaller time step, At = 0.05. Although the beat 
pattern is slightly enhanced, likely indicating that we have started further away from the true oscillon because of the 
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FIG. 10: Energy in a box of radius 28 as in Fig. [TJ but with initial conditions that have been deformed to break rotational 
symmetry. The spatial coordinate r = xx + yy + zz has been everywhere replaced by r' — 0.98xx + 1.02yy + 0.97zz and 
similarly r and r have been replaced by r' = jr'j and r' — x' /r' . As a further check of the numerics, this run also uses a smaller 
time step, At = 0.05. Except for these modifications, the simulation is the same as in Fig. [T] 



nonsphcrical deformation, we see that the system nonetheless converges to a very similar configuration to the case 
without the rescaling. Equivalent behavior is seen when we make these two changes individually and when we make 
other perturbations, such as variations of x- 



Finally, we consider the topological properties of the electroweak oscillon. Unfortunately, as shown in [47j, there is 
no unambiguous definition of the topological charge for solutions to the equations of motion. (Topological properties 
are typically studied using vacuum-to- vacuum paths [48} , which are clearly not solutions to the equations of motion 
since they do not conserve energy.) However, for any localized spatial configuration in which the Higgs field never 
vanishes, the Higgs winding number is unambiguously defined as 

jk Tr [U^diUp^djUp^dkU)] d 3 x , (27) 
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where U is the unique SU(2) matrix associated with a nonvanishing Higgs field so that 

$ = \$\u(iy (28) 

The Higgs winding number is a topological invariant, which can only change with time if the Higgs field passes through 
zero at some point in space. The change in the Higgs winding is physically meaningful and measures whether the fields 
have crossed the sphaleron barrier. Because the electroweak oscillon contains only small-amplitude field fluctuations, 
its Higgs winding is always zero and it does not approach the sphaleron barrier. Correspondingly, its topological 
density 

q = ^ VX ° F ^ ■ (29) 

is small as well. But the restriction to small amplitude does not apply to its decays (induced, for example, by 
collision with another oscillon), when the fields frequently exhibit an implosion to small radii and large amplitudes 
before ultimately dispersing. This behavior is seen in Fig. [3] for the oscillon's decay when A is slightly less than one. 
However, both this particular decay and limited experiments with oscillon collisions have not led to winding in the 
final Higgs field. Current work continues to investigate this possibility. 



CONCLUSIONS 



We have seen in detail the results of a numerical simulation describing a long-lived, localized, oscillatory solution to 
the equations of motion in the bosonic sector of the electroweak Standard Model, for a Higgs mass that is twice the 
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W mass. Compared to the natural scales of the system, this solution has small field amplitudes, large spatial extent, 
and large total energy. In the quantized theory, it would represent a coherent superposition of many elementary 
particles, and thus is well described by the classical analysis undertaken here. Quantization of the small oscillations 
around the classical solution would nonetheless be of interest, as has been done for Q-ball oscillons in [49]. It would 
also be desirable to incorporate fermion couplings, which have been ignored here. Such an analysis would require 
introducing chiral fermions on the lattice, which is well known to be a difficult problem, but one on which significant 
progress has been made in recent years. While one might expect the oscillon to be destabilized by decay to light 
fermions, in the case of the photon coupling we have seen that the analogous decay mechanism is highly suppressed. 

Because it would require bringing many Higgs and gauge particles together at once, forming such an oscillon 
would likely require large energies available only in the early universe. If extremely long-lived, such an oscillon could 
be a dark matter or ultra-high energy cosmic ray candidate. A slow fermion decay mode would be of interest for 
baryogenesis, since it could provide a mechanism for fermions to be produced out of equilibrium, as is necessary to 
avoid washout of particle/antiparticle asymmetry. The oscillon has small amplitude everywhere and thus remains 
far from the sphaleron configuration, even though it has energy above the height of the sphaleron barrier. However, 
when induced to decay, for example by a collision with another oscillon, the fields typically collapse to a configuration 
with small radius and large energy density and field amplitudes before dispersing. Such decays could potentially cross 
the sphaleron barrier and produce fermion number violation. For baryogenesis applications, one would also need to 
incorporate interactions containing C and CP violation in the classical effective action. 

The spherical ansatz provided a crucial tool for obtaining the electroweak oscillon solution. However, any search 
for oscillons using a particular ansatz cannot guarantee that all solutions have been found. "Emergent" techniques, 
in which oscillons form from generic initial conditions, offer the opportunity for more comprehensive searches for 
oscillons, albeit at a higher computational cost. In simpler models, oscillons have been shown to emerge from phase 
transitions [30| and from thermal initial conditions in an expanding universe [Ho| . Clearly, it would be desirable to 
extend these techniques to the electroweak model. 

The electroweak oscillon remains stable even when one would expect it to decay, suggesting that there might exist 
other stable, oscillatory solutions in the electroweak theory or its extensions, either for generic or specific mass ratios. 
While results for generic mass ratios are clearly of broader applicability, a compelling result for a specific mass ratio 
might suggest a preferred value of the Higgs mass. 
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